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> Magnetocardiographic (MCG) 
signal is a result of 
simultaneous electrical activity 
of masses of heart cells 


HBasic morphologies of the 
MCG signal are equal with 
ECG signal 



Signal consist on 3 basic parts; 
P wave, QRS complex and 
T wave 



QRS 
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Recording of the MCG signal 


Biomag2000 


Signal is recorded, usually in 
magnetic shielded room, with 
SQMDs 

Analog signal is low-pass filtered 
and digitized 

Most devices nowadays are 
multichannel devices, enabling 
mapping studies 




Example of measurement 
grid (Biomag laboratory) 
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Preprocessing of the data 


Signal is triggered 
► Usually from R peak 

Baseline is detected and subtracted from signal 
Signal is averaged 
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Fitted spline (3-deg) baseline 
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Analysis based on source localisation 


> Ill posed problem 

► Some restrictions for the source 
are needed 

Point source models (current 
dipoles) 

For localized sources 

M Distributed source models 

For distributed or propagating 
sources (activation propagates as 
wavefronts) 



Typical field and 
potential maps 

















> Multichannel MCG data 
from the beginning of the 
ectopic beat are used for 
localization 

M Torso is modelled with 
boundary elements 

H Position and direction of 
the dipole are iteratively 
evaluated 



Localization of arrhythmic foci 
(20 points) 
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Analysis based on source localisation (current distribution) 

Localisation of the injury current 


m 


> Data are the difference of the 
ST-level before and after 
exercise test 

>j Current density estimate 

m The magnitudes of the dipoles on 
preselected surface are evaluated 
by solving the inverse problem 

H The result is visualized 
(on 3D or bull’s eye view) 


Ant 



Max 

CDE p °s 



Bull’s eye view of the 
injury current estimate 

[K. Pesola 1999] 
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Analysis based on source localisation (current distribution) 

Activation time map 



The positions and the 
orientations of the dipoles 
are preset (uniform double 
layer) 

The sequence of 
depolarization during the 
QRS is solved from 
measured multichannel MCG 
data 




Invasive 

method 



Activation times 

[T. Oostendorp 1998] 
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Analysis based on source localisation 

Localization with QRS integral maps 

Mean QRS integral maps during the induced VT 
and ventricular pacing are evaluated (from BSPM ) 

[B. Horacek 1979] 

Pacing localization is changed until the maps equal 



Mean QRS integral maps over the shaded are 
(left: induced VT, right: ventricular pacing ) 
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Is Analysis based on source localisation 

Classification of maps with SOM and LVQ 

> SOM was trained with potential (field) maps taken 
from different patients at different time instants 
during the QRS 

>1 After training the SOM the distance maps or the 
trajectories are calculated and classified 
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Distance maps on SOM 


L. Reinhardt 1998 



SOM trajectory 

L. Reinhardt 1998 


11 













Biomag2000 


Analysis of the signal averaged data 


> Basic preprosessing 

> Improves the signal to noise 
ratio 

■ > Emphasis cyclic features 

M Averaging is usually done 
somehow selectively 



An example of averaged MCG data 
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Analysis of the signal averaged data 

P wave 

P wave is analysed in a similar manner than the 
QRS (next slides) 

P wave itself is used as trigger instead of the QRS 



Averaged and non averaged P wave 
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Analysis of the signal averaged data (QRS complex) 

Time-domain analysis 

Late-field analysis (simson method) 

> High-pass filtered (Butterworth 40 Hz) [M. Simson 1983] 

LHEnvelope of the filtered signal is formed using the 

Hilbert transform [J. Montonen 1994] 

Parameters: QRS duration, RMS, LAS 

QRSdur = 140 ms 
RMS40 = 392 fT 
LAS1 pT = 62 ms 


Filtered QRS 



Original signal 
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Analysis of the signal averaged data (QRS complex) 

Fragmentation analysis 

Fragmentation analysis [P. Endt 1998] 

> Intra-QRS analysis 
LaBand-pass filtering 37-90 Hz 
(binomial) 

Parameters: Number of extrema (M), 

Fragmentation score (S) 



Original signal Filtered signal 
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Analysis of the signal averaged data (QRS complex) 

Spectral turbulence analysis 


Spectrotemporal analysis [G. Kelen 1992] 

> Intra-QRS analysis 

LB Quantities the variability between the short time (24 ms) 
Fourier spectra within the QRS 

Parameters: QRS duration, ISCM, ISCSD LSCR, 
SE 
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Analysis of the signal averaged data (T wave) 

QT dispersion 


QT-dispersion analysis 

> Quantities the spatial QT time differences 

Parameters: Max(QT) - Min(QT), 

Max(Tend-Tapex) I Min(T end -T apex ), QT_SD, local QTd 


2 pT 
100 ms 


309 ms 


QRS and QT durations 


QT apex isovalue colour map 


17 








Biomag2000 


Analysis of the signal averaged data (T wave) 

Other methods 

QRST-integral 

Relative smoothness (RS) score 

The 90 % duration of the T wave 

Time instants of maximum onset and offset slopes 
of the T wave 

Symmetry parameter (SR) of the T-wave 
The number of the ‘phases’ in the T-wave 




Dynamical analysis 
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Analysis of interbeat variation or changes of signal 
‘Spontanic activity’ 

J Prefers long term recordings 

I « Heart rate variability, QT dynamics, T-wave alternans, 



Interventions, exercise testing 

Data is recorded before, (during) and after the 
intervention 

Especially used for ischemia detection 
ST-depression, ST- and T-maps, QRS envelope, 
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Dynamical analysis 


■F" 

dl 


Heart rate variability 

Time series must be corrected from ectopic beats 

Time domain analysis 

Average RR time and standard deviation 

Frequency domain analysis 

laPowers of different frequency bands: VLF, LF and HF 


\ w 

50 s 

Ml 

Uncorrected RR timeseries 




20 ms 


50 s 


0.1 Hz 


Corrected timeseries 


Power spectra 
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Dynamical analysis 

QT dynamic 


> QT / RR relation 

► Depends on previous RR intervals 

► Time domain analysis of QT timeseries 
Frequency domain analysis of the timeseries 
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QT/ RR plot (24-h Holier ECG) 


RR time (ms) 

QT/ RR plot with colour coded 
preceding heart rate 
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Dynamical analysis 

T wave alternans 


Consistent 2:1 variation of the 
waveform morphology 
>- Depends on heart rate 

Usually studied in frequency 

domain [D. Rosebaum 1994] 

1, The power spectrum of 
consecutive T wave amplitudes 
is formed 

2. The frequency 0,5 is scaled 
with noise value 


CD 


Alternans ratio =1.12 
Cum alt voltage = 0,64 




0,0 0,25 0,5 

Frequency (1/beat) 

Alternans spectrum 
(24-h Hotter ECG) 
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Dynamical analysis (interventions) 

Envelope of the QRS 

Detection of high frequency signal attenuation 

[S. Abboud 1992] 

> The QRS complex is band-pass filtered (150-250 Hz) 

IflThe envelope of the filtered signal is formed 

I ^[Difference between the envelopes before and after the 
•.interventions are calculated 



20 ms 

Filtered signals (before (black) and Difference of the envelopes 
during exercise (white)) 
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Dynamical analysis (interventions) 

Rotation of T and ST integral maps 


1. The integral field maps of either ST -segment or T wave 
are evaluated 

2. The angle of the maximum derivative is determined 
from the integral field maps 

3. The angle is plotted as the function of time or as the 
function of heart rate 
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[P. Takala 2000] 


[P. Takala 2000] 
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> There is no ‘ultimate method’ 

> Different analysis methods for different tasks 

> Many ways to analyse one data 

Method to efficiently combine the results of 
different analyses is needed 




